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DEPOLARIZED  INFRARED  REFLECTANCE 
FROM  DRY  AND  WETTED  SURFACES 


X.  INTRODUCTION 

In  this  study,  we  incorporated  two  experimentally  verified  properties*  of 
light  scattering  to  the  active,  direct  detection  of  a  targeted  liquid  on/within  sin 
inhomogeneous  dielectric  surface.  One  property  is  that  multiple  scattering  is  a 
subsurface  effect  and  accounts  for  depolarized  radiance;  the  other  property  is  that 
this  depolarized  component  is  Lambertian  (isotropic  radiance)  and  not  correlated  to 
topographical  detail  of  the  irradiated  surface.  The  goals  of  these  investigations  are 
to  separate  and  measure  volume -scattering  powers  and  to  determine  the  spectral 
sensitiveness  versus  applied  contamination  density.  Although  volume  scattering  in 
the  mid  infrared  (MIR)  accounts  for  only  a  minor  percentage  of  the  total  scattering 
signal  from  soil,  sand,  asphalt,  and  concrete  surfaces  (these  surfaces  of  interest 
absorb  MIR  radiation,  and  the  stronger  the  absorption  the  less  the  volume 
scattering),  we  suspect  that  this  component  alone  contains  a  greater  percentage  of 
information  influenced  by  the  liquid  and  that  by  being  a  subsurface  effect  is 
insensitive  to  surface  roughness.  The  work  of  Renau  et  ah*  was  extended  from  the 
visible  wavelength  of  0.6328  nm  to  the  9-  to  ll-jim  band  of  the  MIR,  and  the 
dielectric  media  were  contaminated  by  ejection  of  a  nonvolatile  liquid  with  strong 
resonant  absorption  centered  about  9.7  fim.  As  a  consequence  of  the  greater 
wavelengths  of  our  probe  beam,  the  depolarized  component  makes  up  only  about 
10%  of  the  total  backscattered  light  in  comparison  to  about  80%  from  the  scattering 
of  a  He— Ne  beam  by  soil  particles,  for  instance.  The  liquid  lying  on  and/or  diffused 
through  the  surface  is  targeted  by  its  ability  to  absorb  in  the  MIR  band.  Selective 
absorption  is  the  basis  for  discriminating  against  the  target  by  irradiating  the 
contaminated  surface  at  alternate  wavelengths  centered  on  and  off  the  absorption 


peak  of  the  liquid  (laser  line  bandwidth  <<  liquid  absorption  bandwidth)  and 
analyzing  the  difference  in  backscattered  signals.  After  filtering  out  the  reflectance 
spectral  response  of  the  unwetted  surface,  a  condition  can  be  defined  that  will  cause 
an  alarm  to  trip,  for  instance,  based  on  the  magnitude  in  the  difference  signals, 
indicating  the  probability  of  an  absorption  match  between  sample  and  target. 

2.  THEORY 

2.1  Liquid  Coating  Formed  on  an  Impermeable  Volume  Scattering 

Substrate. 

A  restrictive  class  of  contaminated  surfaces  are  those  that  will  retain 
dispersed  liquid  droplets  on  their  surface.  A  metal  plate  with  military— specified 
primer  coating  and  paint  layer,  which  was  used  in  these  experiments,  is  such  a 
surface.  After  striking  the  paint  layer,  the  droplets  spread  and  eventually  form  a 
thin  film.  The  method  for  detecting  this  liquid  layer  is  shown  in  Figure  1.  A 
linearly  polarized  incident  beam  generates  a  source  of  volume -backscattered 
radiance  that  back— transmits  through  the  liquid  layer.  The  volume  reflectance 
signature  will  therefore  maximize/minimize  at  wavelengths  that  axe 
nonabsorbing/absorbing  by  the  liquid.  This  differential  reflectance  feature  contains 
quantitative  (possible  density  mapping  applications)  as  well  as  qualitative 
(detection)  information,  since  the  depolarized  reflectance  component  is  related  to 
film  thickness  (Equation  3b)  via  amplitude  attenuation  of  the  multiple— scattered 
waves  transmitting  through  the  liquid  layer  at  wavelengths  where  refractive  index 
in  the  liquid  is  complex  (k  ^0). 


Depolarized  Reflectance  Modeling  Equations. 


From  specular  reflection  properties  of  the  liquid  layer  and  radiometric 

3 

definitions  of  radiance  and  irradiance,  a  volume  reflectance  function  that  solves  the 


FIGURE  1 .  The  proposed  concept  for  djrect  detection  of IR  absorbing  liquid 
on  an  impermeable  dielectric  surface.  k;n,  ksp,  &ss  and  kms  are  wave  vectors 
of  the  incident  parallel-polarized  wave;  the  specularly  reflected  parallel- 
polarized  wave;  a  modified  single-scattered  parallel-polarized  wave;  and  a 
modified  multiple-scattered  polarization-shifted  wave,  respectively.  P  is  a 
polarizer  with  transmission  axis  perpendicular  to  the  plane  of  incidence. 
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I*. 


boundary  value  problem  depicted  in  Figure  1  was  obtained  and  is  given  below.  This 
solution  is  a  variation  of  the  model  given  in  Wolfe  and  Zissis’  handbook^  rederived 
by  Rosengreen  and  coworkers'*  and  is  nonrigorous  (pv  is  an  experimental  input) 
with  the  assumption  of  a  priori  Lambertian  volume  reflectance  from  the  unwetted 
surface.  The  model,  however  limited  in  usefulness,  can  provide  relative  reflectance 
differences  as  modified  by  varying  layers  of  the  absorbing  film.  The  selection  of 
wavelengths  that  most  discriminate  against  the  contaminant  is  thereby  simplified. 
Volume  reflectance  defined  as  the  ratio  of  depolarized  radiant  power  transmitted 
through  the  liquid  layer  to  total  incident  beam  power  (Figure  1)  is  given  by: 

fv=Pv - 7— —  T~~~  (!) 

1-*>*(JnkTdtlh)  *•/* 

where ; 

pv  denotes  volume  reflectance  of  the  unwetted  sample , 
pt  denotes  specular  reflection  of  the  incident  beam  at  the  air —to— liquid  interface, 

*1 

plane  of  incidence.  (2) 

2  II 

t  and  r  are  transmission  functions  given  by : 


cos29i—(n2+k2)cos29t)2+4k2cos29tcos29i 

[( neos  0, +0080^ ) 2 + k  2coe20,- ] 2 

[(fc2+n2)  cos20,-cos20(]2+4k2cos20<cos20, 
[(neos  9i+cos9t)2+k2cos29i)2 


ere 


^5* 


T(ffi)  =  « 


CO80; 


(3a 


r(0,)  =  c 


cos©,- 


(3b) 


and 


_  _  4nwk[X) 
P=  X  ’ 


p,  denotes  reflectance  of  depolarized  radiance  at  the  liquid —to— air  interface, 


6,,  Ot  are  incident ,  transmitted  beam  angles  at  the  air  —to— liquid  interface  , 
©,•,  and  0C  are  incident,  transmitted  and  critical  ray  angles 


at  the  liquid— to— air  interface. 


(  all  angles  are  measured  relative  the  liquid  surface  normal  ) 

and  E,  is  the  Exponential — Integral  function.  ' 

A  wavefront  propagating  within  the  liquid  medium  at  wavelengths  where  k  ^  0  is 

‘inhomogeneous’  by  nature,  viz,  its  surfaces  of  constant  phase  and  constant 

amplitude  do  not  generally  coincide.  This  implies  a  complex  interpretation  of 

Snell’s  law  and  further  implies  that  the  propagation  direction  of  a  plane  wave 

7 

refracted  into  the  absorbing  liquid  layer  is  given  by: 


y  smiOi-\-qi(ncos^— ksin'f) 


where ; 


g2(0,)  =  y/(l  -Asin  20j )  2+(Bsin  i0i ) 2, 


n2[l  +  (k/n) 


n2- e 


(n2+fc2)2  ’ 


*  Internal  3.351  applies  to  Equation  4a  when  substituting  y  =  (co»OJ  .  See  section  8.2  for  properties  of  the 
exponential-integral  function. 


BsinHi 
1 —Asm2#, 


(5e) 


li»i)  =  |tan' 


and  like  expressions  for  a  polarization— shifted,  volume— scattered  ray  inside  the 
liquid  propagating  toward  and  refracted  at  the  liquid— to— air  interface.  Note  that 
Equation  5a  reduces  to  the  more  recognizable  sin#,-  =  nsin6t  Snell  law  when  k=0. 

The  perpendicular  and  parallel  components  of  specular  reflectivity  at  the 
liquid— to— air  interface,  pper  and  ppor  in  Equation  4a,  are  obtained  by  replacing 
(6j,6t)  with  (0(,©,)  in  Equation  2,  and  the  critical  angle  ©c  is  determined  through 
substitution  of  (0C,0(  =  tt/2)  for  (#  t  ,0, )  in  Equations  5a— e.  To  evaluate  Equation 
4a,  a  substitution  of  variable  via  Snells’  transformation: 
0,  -»  g{Qt,n,k);  dsin 0,  — ►  h(Qt,n, k)dQt,  was  required.  The  resulting  lengthy 
expression  for  p„  is  numerically  integrated  by  a  Simpson’s  rule  subroutine  in  the 
model  code. 

The  (1— p,)t  term  from  Equation  1  is  a  measure  of  how  much  energy  is 
available  for  volume  scattering,  while  the  P„{\.-pt)T  term  gauges  how  much  of  that 
energy  transmits  the  liquid  layer  and  is  subsequently  available  for  detection.  The 
denominator  of  Equation  1  results  from  an  infinite  number  of  internal  reflections 
within  the  smooth,  continuous  liquid  layer,  and  by  applying  the  binomial  expansion: 

(l-x)~l  =Ji',  where  x  =  pv{j  Tdnh)2~pjn  «  1. 
i=o 

All  the  ingredients  for  evaluating  Equation  1  from  input  parameters 
A,  n(A)  and  fc(A),  p„(A),  0,-,  Qt,  and  liquid  thickness  w  are  now  provided.  The 
modeling  software  produces  three  dimensional  graphics  for  generating  surfaces  of 
/„(A,u;).  An  /„(A,u;  =  u;0)  cross  section  of  the  computer  generated  surface  can  be 
compared  to  experimental  data  from  the  contaminated  surface.  Detection  is 


L.»  <.«  I.-J, 


».•  I  ,* »  jjjjWjJj  HU  »JL ^  tf .  ^ 


probable,  and  film  thickness  (u>0)  is  extracted  provided  a  spectral  match  condition 
can  be  established. 

2.2  Liquid  Diffused  Into  the  Bulk  of  a  Porous  Volume  Scatttrer. 

A  nonvolatile  aerosol  simulant  (i.e.,  a  liquid  with  rheological  and  optical 
properties  that  approximate  a  toxic  agent  compound  and  is  nonchemically  reactive 
with  the  scattering  particles)  ejected  onto  terrains  will  generally  diffuse  into  a  solid 
bulk  with  little  or  no  surface  retention.  The  interstitial  structure  of  liquid— coated 
particles  and  voids  constitute  a  chemically  complex  scattering  medium  to  which  a 
true,  tested,  rigorous  statistical  reflectance  theory  does  not  exist  to  our  knowledge. 
We  do  not  offer  solutions  to  this  problem,  but  do  suggest  that  the  effective  medium, 
extended  effective  medium,  and  multiple  scattering  theories  of  Bohren  and 
Huffman,  Chylek  and  Srivastava,  Stroud  and  Pan,  and  Veradan  et  ah  are 
useful  for  future  modeling  work.  Before  invoking  the  effective  medium  theories, 
however,  one  must  be  aware  of  their  constraints,  assumptions,  and  limitations.*1 
We  finally  note  that  work  by  Bahar,1^’1^’1^’1^  using  a  ‘Unified  Full  Wave’ 
approach  for  solving  problems  of  scattering  by  nonuniform  stratified  media,  is 
currently  ongoing  under  contract  with  CRDEC,  and  this  may  provide  insight  into 
the  problem  of  the  permeable  contaminated  surfaces. 

3.  EXPERIMENTAL  APPROACH 

The  optical  arrangement,  feedback  control,  data  acquisition  and  processing 

* 

electronics  of  the  experiment  are  schematically  drawn  in  Figure  2.  SF96,  a 
nontoxic  silicon— based  oil,  with  double— peaked  resonant  absorption  profile  between 
9  and  11  ,um,  Figure  3,  served  as  the  contaminant  and  target.  A  COg  waveguide 

*  SF96  is  a  simulant  of  chemical  agent  VX  widely  used  in  military  research. 
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cw  laser,  grating  tunable  from  its  two  P—  and  R-branch  vibrational-rotational 

transitions,  provides  a  selection  of  75  wavelengths  between  9-  and  11— /xm, 

emissions  within  the  MIR  window  of  atmospheric  transmission.  The  CC^  beam  is 

polarized  parallel  to  the  sample  plane,  with  an  extinction  ratio  better  than  100:1, 

and  expanded  to  produce  an  elliptical  cross-sectional  area  of  6.2  sq  cm  for  a  fixed 

50  degree  angle  of  incidence.  Speckle  interference  is  integrated  through  beam 

expansion  and  by  rotating  the  sample  at  1  Hz.  Scattered  radiance  is  directed  to  the 

detector  through  a  three— mirrored  goniometer  arm  that  scams  the  polar  angle  in  the 

plane  of  incidence  across  the  upper  hemisphere.  The  light— receiving  system 

consists  of  a  liquid— nitrogen— cooled,  1  mm  x  1  mm  HgCdTe  photoconductive 

detector  with  an  AC  coupled  preamplifier  of  375  KHz  bandwidth  and  a  detectivity 

of  1.8x10^  cmHz^^/watt,  an  instantaneous  field  of  view  of  0.063  rad,  and 

numerical  aperture  of  1.59.  Incident  beam  power  was  fixed  at  nearly  1  watt  after 

tuning  between  laser  lines  (via  variable  ND  filter)  so  as  not  to  exceed  the  detector’s 

dynamic  range.  Given  these  system  parameters  and  a  typical  volume  reflectance  of 
—3 

10  ,  the  maximum  range  for  noise -limited  signal  detection  was  computed  to  be 

16 

100  meters.  Polarization  states  of  scattered  light  were  ascertained  by  insertion  of 
a  quarter— wave  plate  before  the  polarizer  and  appropriate  adjustment  of  each 
optic’s  axes.  Contaminant  and  target  SF96  was  quantitatively  dispersed  to  the 
test  surfaces  through  am  ultrasonic  nebulizing  nozzle  contained  inside  a  pneumatic 
chamber.  Software  was  written  to  control  rotation  of  the  Ge  polarizer’s  optical  axis 
and  extract  minimum  and  maximum  detector  voltages  for  each  wavelength  in 
measurement  sequence  of  reference,  dry,  and  wet  samples.  The  depolarized 
reflectance  component  is  proportional  to  the  ratio  of  this  minimum  signal  of 
scattered  radiance  to  the  voltage  generated  by  the  incident  split  beam  striking  a 
pyroelectric  detector.  These  ratioed  voltages  were  further  ratioed  to  reflectance 
from  a  standard  Lambertian  reflector,  with  properties  of  spectral  flatness  between  9 


and  11  fim  and  diffuse  reflectance  comparable  to  the  sample.  (‘Flowers  of  sulfur’ 
powder  is  a  common  standard  used  in  IR  spectroscopy.)  Absolute  reflectance  is 
obtained  from  a  separate  experiment  by  referencing  the  standard  reflector  to  a 
known  gold  reflector. 

4.  RESULTS,  DISCUSSION,  CONCLUSION,  AND  FUTURE 
WORK 

Scattering  by  all  but  one  surface  produced  a  combination  of  linear  and 
depolarized  components.  The  exception  was  a  laboratory— grade  powdered  sulfur 
sample  that  produced  depolarized  and  elliptically  polarized  components  because  of 
birefringence  in  the  sulfur  crystals.  This  powder  was  originally  intended,  but  later 
rejected,  for  use  as  a  standard  reflector.  The  polar  angle  dependence  of  volume 
scattering  satisfied  the  cos<?  Lambertian  condition  from  all  measured  dry  surfaces 
without  exception.  By  experiment,  we  attempted  to  establish  the  relationship  of 
volume  reflectance  with  the  ratio  of  particle  size  of  a  dielectric  scatterer  to 
wavelength  of  the  incident  beam.  The  experimental  results  have  shown  that  the 
depolarized  power  component  increased  by  a  factor  of  ~3  between  dry  alumina 
crystals,  with  particle  size  distributions  centered  at  A  and  A/5,  and  by  a  factor  of 
«40  between  A  and  A/25  crystals.  In  general,  depolarized  power  backscattered  from 
ail  irradiated  dry  materials  prepared  for  these  experiments  exhibited  spectral 
dependencies  -  A  n.  Paint  and  soil  samples  (Figures  4a  and  5a,  dry  data),  for 
instance,  have  shown  Rayleigh -like  scattering  behavior,  where  n  was  calculated  by 
least—  squares  analyses  to  be  4.3  and  4.5,  respectively. 

Distinct  MIR  spectral  differences  between  dry  and  SF96-sprayed 

impermeable  paint  samples  were  observed  to  densities  of  the  contaminant  as  low  as 
2 

2  gm/M  ,  as  illustrated  by  the  typical  measurement  graphed  in  Figures  4a, b.  From 
the  limited  reflectance  data  sets  measured  from  dry  and  wetted  soil,  sand,  concrete, 
and  asphalt  samples,  we  generally  conclude  that  the  spectral  response  of  the 
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FIGURE  4.  Volume  reflectance  (a)  and  corresponding  difference  (b)  spectra 
between  a  dry  (O)  and  SF96-wetted  metallic  slab  with  IR  primer  and  paint 
layers.  The  densities  of  dispersed  contamination  are:  10(d).  5  (  o  ).  2(a) 
gm/M2.  An  analytical  (9P28, 9.621  /<)  and  reference  (10P18.10.571  u)  set  of 
probe  emissions  are  indicated  in  Figure  b. 
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FIGURE  5.  Volume  reflectance  (a)  and  corresponding  difference  (b)  spectra 
between  dry  ( O )  and  SF96-wetted  soil  samples  contaminated  to  densities  of 
50(0)  and  100  (V)  gm/M2.  An  analytical  (10P18.10.571//)  and  reference 
(9P28, 9.621//)  set  of  probe  emissions  are  indicated  in  Figure  b. 


depolarized  reflectance  component  from  permeable  non  homogeneous  dielectric 
materials  is  less  sensitive  to  the  SF96  target  by  at  least  an  order  of  magnitude  of 
liquid  mass.  As  an  example,  our  results  for  dry  and  contaminated  soil  are  shown  in 
Figures  5a, b. 

The  /„  model  code  developed  for  the  liquid  -  coated  surface,  Figure  1,  was 

applied  and  compared  to  the  experimental  results  in  Figure  4a.  The  3  dimensional 

graphics  output  was  successful  at  predicting  the  extremum  located  at  10.57  pm  and 

how  thick  an  SF96  film  was  formed  after  aerosol  ejection  onto  the  painted  metal 

surface.  MIR  absorption  through  the  film  layer  is  apparent  from  these  data  to  a 

2 

contamination  density  of  2  gm/M  ,  perhaps  slightly  less.  Although  this  may  seem 

encouraging,  no  claim  is  made  that  the  data  can  uniquely  determine  the 

contamination;  namely,  the  probability  of  false  alarming  can  be  minimized  or 

eliminated.  Indeed,  it  is  not  unique  since  nonvolatile  liquids  with  overlapping 

absorption  bandheads  (  active  detection  on/within  surfaces  must  successfully  target 

a  class  of  organophosphonate  contaminants,  high  molecular  weight  liquids  that 

generally  exhibit  broad  MIR  absorption  bands)  can  produce  inseparable  differential 

reflectance  signatures.  As  previously  mentioned,  the  contaminated  porous  surfaces 

are  theoretically  difficult  to  model  and  experimentally  much  less  spectrally  sensitive 

to  the  diffused  contaminant.  A  constraint  imposed  on  sensing  the  diffused  liquid 

within  the  bulk  of  impermeable  terrains  using  differential  volume  scattering  is  the 

depth  to  which  the  probe  beam  can  penetrate  the  surface  and  thereby  interact  with 

the  liquid  cells.  Recall,  this  is  one  reason  volume  scattering  was  separated  and 

I ')  * 

chosen  as  a  feature  for  discrimination  From  optical  constants  data  of 
montmorillonite,  kaolinite  and  illite  compressed  slabs,  basic  minerals  composing  soil, 
a  penetration  depth  at  the  e  ^  point  minimizes,  unfortunately,  to  alxuit  I  pm. 
between  wavelengths  of  9  and  10  pin 
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Mid  infrared  volume  reflectance  from  dry  inhomogeneous  dielectric 

materials  was  found  to  be  isotropic  and  inversely  dependent  on  wavelength  to  some 

power.  This  depolarization  component’s  sensitivity  to  an  MIR-absorbing  liquid 

diffused  through  porous  surfaces  is  limited;  the  lower  limit  %  20  gm/M  SF96 

sprayed  onto  terrains.  For  an  impermeable  surface  boundary,  however,  the  data 

o 

suggests  that  SF96  absorbency  can  be  revealed  at  densities  as  low  as  2  gm/M  (a  2 
Jim  film).  These  findings  may  spur  continuing  efforts  on  direct  detection  methods 
to  improve  sensitivity,  however,  new  technologies  which  are  more  sensitive  to  and 
more  fully  characterize  the  scattering  surface  must  be  explored  to  fulfill  the  Army’s 
remote  sensing  mission.  Current  and  future  experiments  in  this  and  other 
laboratories  involve  exploration  of  other  light  scattering  technologies  that  may 
provide  greater  sensitivity  and  precision  (some  toxic  liquids  of  interest  and  their 
non  toxic  derivatives  do  exhibit  overlapping  absorption  bandheads  between  8  12 
jim)  in  discriminating  against  surfaces  targeted  for  known  hazardous  contaminants. 
Current  work  by  one  of  us  (Carrieri)  is  directed  toward  measurement  and  analysis 
of  all  16  elements  of  the  Mueller  matrix,  the  4x4  transformation  between  incident 
and  backscattered  Stokes  vectors  that  fully  characterizes  the  scattering  surface  as 
functions  of  scattering  angle,  contamination  density,  and  wavelength  between  the 
9  12.5  jim  ^C^O^,  12CI8Q2,  and  ^C^Ot>  laser  emission  bandwidth. 
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